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e Abstract 

7 We report experiments on the production of ultracold neutrons (UCN) in a converter of 

8 superfluid helium coated with fluorinated grease. We employed our technique of window-free 

9 extraction of accumulated UCN from the helium, in which they were produced by downscat- 

10 tering neutrons of a cold beam from the Munich research reactor. The time constant for UCN 
n passage through the same extraction aperture as in a previous experiment was a factor two 

as : 12 shorter, despite a lower mean velocity of the accumulated UCN in the present experiments. 

13 A time-of-flight measurement of the cold neutron spectrum incident on the converter allowed 

14 us to estimate the multi-phonon contribution to the UCN production. The UCN production 

15 rate inferred from two methods agrees with the theoretical expectation. 

o 

16 PACS numbers: 78.70.Nx, 28.20.Fc, 29.25.Dz, 61.12.Ha 

17 Keywords: ultracold neutrons, UCN, UCN sources 



*email: zimmer@ill.fr 



1 Introduction 



20 Ultracold neutrons (UCN) have energies in the neV range and velocities up to a few meters per 

21 second. When impinging on suitable materials they undergo total reflection under any angle 

22 of incidence and can therefore be trapped in bottles and manipulated for a long time (see the 

23 books [1, 2] for an introduction to the physics of UCN). Owing to this feature UCN have become 

24 very useful in various fundamental investigations of neutron properties, with strong implications 

25 for particle physics and cosmology. The longstanding search for the neutron electric dipole 

26 moment investigates CP-violation beyond the standard model of particle physics [31 0]. The 

27 accurate determination of the neutron lifetime is required for a detailed understanding of big 

28 bang nucleosynthesis [5]j an d to investigate strength and structure of the semi-leptonic weak 

29 interaction within the first quark family (see, e.g., the workshop proceedings [6j [T] ) . New topics 

30 with UCN are, among others, the demonstration of quantum levels of neutrons in the earth's 

31 gravitational field [8 J , and even more recently, the search for neutron - mirror neutron transitions 

32 [91 [10] . The most intense UCN source at the Institut Laue Langevin in Grenoble [TTj provides 
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33 not more than 50 UCN per cm 3 . In order to improve counting statistics, new UCN sources are 

34 being developed in many laboratories around the world [I2j [I3j [HI [151 [I6j fT7[ [18] . 

35 An elegant method to produce UCN employs a converter of superfluid 4 He in a beam of cold 

36 neutrons [19J. The kinematics defined by the dispersion relations of helium and the free neutron 

37 enables downscattering of cold neutrons with an energy around 1.0 meV (wavelength 0.89 nm) 

38 to ultracold energies via emission of a single phonon. In addition, multi-phonon processes occur 

39 which contribute to the integral UCN production rate for a wide range of incident neutron 
« energies [201 [21]. At low temperatures the probability for upscattering is strongly suppressed by 

41 the Boltzmann factor. Therefore, and since pure 4 He has no cross section for neutron absorption, 

42 a large density of UCN may build up in a converter with reflective walls. The storage time 

43 constant r can of course not exceed a limit close to 900 s set by the neutron beta decay lifetime. 

44 In a recent experiment we have demonstrated for the first time that one may efficiently extract 

45 UCN from the converter after having them accumulated therein [22J. As in past experiments 

46 [161 [T71 [23l [24"] the UCN production rate was found to agree reasonably well with the theoretical 

47 expectation. We have developed a new cryostat designed to keep the source portable and easy 

48 to operate. With a short cooling cycle of a few days, our system is much more flexible than 

49 an earlier apparatus [25] . which was designed to be installed close to the target of a spallation 
so source. Our present apparatus is a prototype for a future UCN source to be installed at a strong 

51 cold neutron beam of a high-flux reactor, where no extraordinary cooling power is required. 

52 A particular feature is the vertical extraction of the UCN through a cold mechanical valve 

53 situated above the helium bath. In contrast to a previous attempt to extract accumulated UCN 

54 horizontally [26], no gaps or windows are required in our method. To gain first experience with 

55 this system, the UCN converter vessel and extraction guide system were made of stainless steel. 

56 In the first run it enabled us to measure, with negligible background, the UCN production rate 

57 and to study the temperature-dependent storage properties of the converter. Here we report 

58 new results obtained with the converter vessel coated with fluorinated grease (Fomblin). It has 

59 good reflection properties for UCN and, to our very surprise, positively influences the extraction 

60 time constant. Moreover, a time-of-flight (TOF) measurement of the incident cold neutron beam 

61 helped to improve the comparison with the theoretical UCN production rate. 

62 2 Apparatus 

63 The apparatus was already described in some detail in ref. [22]. The converter vessel has a 

64 volume of about 2.4 liters, made from an electropolished stainless steel tube with length 696 

65 mm and inner diameter 66 mm, closed by Ni windows on both ends. The lowest temperature 

66 attained for the completely filled vessel was T = 0.72 K in the previous run, and 0.82 K in the 

67 present. The maximum kinetic energy of UCN storable in the vessel is defined by the Fermi 

68 potential Vp jWa n of the wall material, which is (184 ±4) neV for the stainless steel used (and 252 

69 neV for the Ni windows), from which one has to subtract the Fermi potential of the superfluid 

70 helium (V F 4 He = 18.5 neV). Due to a small neutron absorbing aperture with diameter 33 mm 

71 placed at the entrance window to the converter, chosen in order to avoid activation of the 

72 vessel in these first test experiments, the volume intersected by the cold neutron beam ("UCN 

73 production volume") was V p = 595 cm 3 , only In the experiments reported here, the inner 

74 surface of the vessel was coated with a thick layer of fluorinated grease (Fomblin) with a Fermi 

75 potential of (115 ± 10) neV. UCN were extracted as previously through a flapper valve situated 

76 above the superfluid helium in the "T" section of the storage tube, connecting the (uncoated) 

77 extraction line made of electropolished stainless steel to a 3 He-gas UCN detector. The present 

78 experiments were again performed at the neutron guide "NL1" at the Munich research reactor 

79 FRM II, using the same neutron beam collimation as in the previous setup [22] . 
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Figure 1: Time histograms of UCN count rates, measured in "buildup mode" (see text) at 0.82 
K for 23 different accumulation times to. 

so 3 Definition of time constants and measurements 

si UCN storage and extraction from the converter can be characterised by various time constants. 

82 The storage time constant r quantifies the temporal decrease of UCN in the closed vessel. The 

83 rate r _1 contains a contribution rj -1 due to wall collisions, absorbing impurities, and UCN 

84 escaping through small holes in the vessel. This contribution does not depend on temperature 

85 T but on the kinetic energy E of the UCN. In addition there are the rates for the T-dependent 

86 UCN upscattering, and for neutron beta decaj0, 

r- 1 =r - 1 ( J B) + r- p 1 (T)+^- 1 . (1) 

87 The emptying time constant r e quantifies the temporal decrease of UCN in the vessel with 

88 the UCN valve open. Therefore, 

t- 1 (T,E)=t- 1 (T,E)+tX 1 (E), (2) 

89 where ta is the time constant for UCN passage through the extraction hole with area A. From 

90 these time constants we may derive the detection probability W for a UCN created in the 

91 converter vessel. It characterises the efficiency of the whole system including extraction and is 

92 given by 

W = e^=e^^, (3) 

93 where the factor e describes losses in the extraction line and imperfect detector efficiency. The 

94 conversion process employed in the present experiments produces a broad spectrum of low-energy 

95 neutrons. Only neutrons with energies below T^ )Wa u — V F 4 He are trapped inside the converter 

96 vessel, whereas those with higher energies quickly escape. The spectrum of the neutrons re- 

97 maining in the vessel is shaped due to wall losses. The energy dependence of these losses is 



1 In the present experiments this contribution forgotten to mention in ref. [22] was still small compared to the 
sum of the other rates. 
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Figure 2: UCN count rate integrals as a function of the accumulation time to- The fit to the 
data points at T = 1.26 K is performed with the single-exponential function defined in eq.Q. 
The results for the fitting parameters r and N are employed in section 4 to determine the UCN 
production rate. The fits to the data points for lower temperatures employ double exponentials 
and serve as guides to the eye. 



98 due to the increase with UCN energy of both the average loss probability per wall collision and 

99 the frequency of wall collisions. Thus the largest losses occur for those neutrons with energies 

100 close to the Fermi potential. However, if on a given time scale the change of the energy spec- 

101 trum stays sufficiently small, the time constant t is still well defined. This is found to be a 

102 reasonably good approximation for the accumulation of UCN at high temperatures, where the 

103 energy-independent (T) dominates the rate r _1 (see eq.flTJ)). At lower temperatures we may 

104 still deduce values for the various rs (still calling them "time constants"), and study how they 

105 depend on the times of UCN accumulation and trapping. 

106 The time constants r and T e can be obtained from measurements in the "buildup mode". 

107 There, the closed converter is first irradiated with cold neutrons for a time to, after which the 

108 beam is shut off and simultaneously the UCN valve is opened. During the whole process a time 

109 histogram of the UCN count rate is recorded. Figure 1 shows a series of such histograms for 
no different accumulation times to, which was obtained for the lowest temperature of the converter 
in in the present runs. The decrease of the UCN count rate while emptying the vessel proceeds with 
ii2 a pure single-exponential; fits of the whole decay in each of the histograms with Aexp(— t/r e ) 
in always resulted in a reduced \ 2 close to unity and provided a value for T e for each of the 

114 histograms. 

115 Figure 2 shows the integrals N (to) of all histograms measured. It demonstrates the satura- 

116 tion behaviour and the strong T dependence of UCN accumulation for the range of temperatures 

117 investigated. If the time constant r is sufficiently well defined, one can fit to the data the single- 

118 exponential buildup function 

N(t ) = N(l-exp(-toM). (4) 

119 For the highest temperature, T = 1.26 K, this fit has a reduced x 2 °f 3.5. The fit becomes 
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Figure 3: Emptying time constants t c as a function of the UCN accumulation time to, determined 
from the histograms of all buildup- mode measurements for the various converter temperatures. 



120 increasingly worse for the lower temperatures. As visible in fig. 3, there is also a dependence of 

121 r e on to j which becomes increasingly pronounced for lower temperatures. These facts indicate 

122 the increasing influence of the first term in eq.JI]) with decreasing temperature. The fit of a 

123 constant value to the data for r e at 1.26 K has x 2 = 2.3. For later use we summarise below the 

124 results of the fits to the data at this temperature, and the corresponding value for ta, obtained 

125 with eq.([2|) and the value for r from eq.©. The uncertainties stated for r, N, and r e have been 

126 scaled to provide a fit with x 2 = 1: 



t = (15.79 ± 0.82) s, 

N = 2740 ±73, (5) 
r e = (9.92 ±0.15) s, 
TA = (26.7 ± 2.6) s. 

127 From the buildup-mode measurements as described before we cannot extract values for r 

128 as a function of the time the UCN stay trapped in the vessel. In order to demonstrate this 

129 dependence we employed a variant of such measurements with a delayed extraction. There 

130 again the converter is first irradiated with cold neutrons while the UCN valve stays closed. 

131 After an accumulation time to the beam is shut, but the UCN valve is opened only after a 

132 delay time t^- A series of such measurements provides histograms for the same to but various 

133 td- A plot of the count rate integrals of two series, for to = 25 s and 250 s, is shown in fig. 4. 

134 Values for r can be obtained from single-exponential fits to a group of several data points. The 

135 variation of slope of the curves shown in fig. 4 indicates a more pronounced increase of r with 

136 £(j when UCN are accumulated for a longer time. This is to be expected, since after to = 250 s 

137 the relative abundance of slower UCN with respect to the (more abundant) faster ones will be 

138 significantly higher than after 25 s, due to the longer storage time constant of the slower UCN. 

139 Their presence in the vessel then becomes better visible in the delayed-extraction experiments. 
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Figure 4: Count rate integrals of the measurements at 0.82 K with delayed extraction after UCN 
accumulation times to = 25 s (lower data points) and 250 s (upper data points). The solid lines are 
fits of single exponentials to the four first, respectively, four last data points. The corresponding 
time constants r derived from these fits are also shown. 



T[K] 


1.26 


1.15 


1.05 


0.93 


0.82 


N c [s- 1 ] 


176 ± 2.6 


227.6 ±2.7 


257.5 ± 1.4 


278.3 ± 2.3 ' 


286.6 ±1.6 


W/e [%} 


37.1 ±3.4 


48.0 


54.3 


58.7 


60.4 



Table 1: Measured T-dependent count rates obtained in continuous-mode measurements, and 
values deduced for W/e (see text). 

140 Values for r can also be deduced from the integral UCN counts for only two different delay 

141 times, using the relation 

^d 2 ~ ^di 

]nN(t dl )-lnN(t d2 y K) 

142 Figure 5 shows such an analysis for the series with to = 250 s. The strong variation of time 

143 constants from one to two minutes demonstrates indeed how much at 0.82 K the spectrum of the 

144 UCN remaining in the vessel is shaped by the energy dependence of the wall collisions during 

145 trapping. 

146 In a third type of measurements, called "continuous mode", we measured the steady state 

147 count rate N c with the converter irradiated for a long time with cold neutrons while the 

148 UCN valve stayed open. The count rates measured are listed in Table 1. In the next sec- 

149 tion we will use N c and W/e at 1.26 K as input to determine the UCN production rate. For 

150 the purpose of illustration we present in Table 1 also values for W/e derived from W{T{) = 

151 PU(1.26K)A' c (Tj)/A r c (1.26K), which assumes that the same UCN spectrum prevails for the dif- 

152 ferent temperatures. Remember that, as obvious from fig. 3, this assumption is only valid as 

153 an approximation. The values stated for the lower temperatures might represent the true value 

154 with a relative error in the order of 10 %. 
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Figure 5: UCN storage time constants deduced from eq. (|6]) , using adjacent data points from the 
measurements with delayed extraction shown in fig. 4. The abscissa is given by the mean of the 
two adjacent delay times, (td«+i + £d;)/2. 

155 4 UCN production rate 

156 The UCN production rate can be inferred from two methods, for which we employ the data at 

157 the highest temperature, T = 1.26 K, where the time constants are sufficiently well defined. 

158 First, as in the analysis of our last experiments [22], we may use the stationary "continuous 

159 mode" count rate N c . The corresponding UCN production rate density Pi is given by 

" * (7) 



V P W 



160 Using the values for N c and W/e stated in Table 1 and the value V p = 595 cm 3 for the UCN 

161 production volume, we obtain 

ePi = (0.797 ± 0.074) s" 1 ™ -3 . (8) 

162 Lacking the knowledge of e, the numerical value in eq.© provides a lower limit for P\. 

163 A second value for the production rate density, P 2 , may be derived from the saturated UCN 

164 number N in the buildup-mode measurement, employing eq.([2J). When UCN are accumulated 

165 for a long time, to 3> T > the production rate equals the UCN loss rate. Denoting with Nq the 

166 saturated number of UCN in the vessel, of which only the fraction N = NqW is detected, this 

167 leads us to 

N 

P 2 = . (9) 

2 WV p t v ; 

168 Using the values for r and N from eq.© with W/e from Table 1 we obtain 

eP 2 = (0.786 ± 0.085) s^cm" 3 . (10) 

169 Hence, within the 10 % accuracy of these two methods, P\ = P 2 - We may now compare the 

170 experimental findings with the theoretical expectation. This requires knowledge of the incident 
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171 cold neutron spectrum which we measured in a separate experiment, using a time-of- flight (TOF) 

172 analysing device described in ref. [27] . It consists of a mechanical chopper and a 3 He detector 

173 with a vertical slit. During a measurement the detector was moved horizontally in order to 

174 integrate over the whole divergence of the beam. The entrance aperture of the TOF analyser 

175 was placed at the position where previously was situated the entrance window to the converter 

176 vessel. We thus determined the neutron wavelength spectrum up to A = 2 nm (frame overlap of 

177 the chopped bunches occured only for A > 5 nm, where the intensity is negligibly small). The 

178 part for A < 1 nm is shown in fig. 6. After the TOF measurements the spectrum was calibrated 

179 with a gold foil activation by the integral neutron flux. 

ibo The mechanism of UCN production contains contributions from single phonon emission and 

181 multiphonon processes. The single-phonon contribution to the production rate density in a 

182 helium converter with Be wall coating is Pi = (4.55 ±0.25) x d0/dAU» s _1 cm -3 , where 

183 the differential flux at A* = 0.89 nm is given in cm _2 s^ 1 nm _1 [17] . For the present situation 

184 this value needs to be corrected for the Fermi potential of the Fomblin grease, i.e. divided by 

185 (252 — 18. 5) 3 / 2 / (115 — 18.5) 3//2 = 3.76, with an uncertainty of 0.59 due to the poor knowledge 

186 of the Fermi potential of the Fomblin grease. From the measured TOF spectrum shown in fig. 6a 

187 we find d(p/d\\ x , = 5.0 x 10 8 cm _2 s _1 nm _1 , from which we expect a single-phonon production 
las rate density Pi = (0.61 ± 0.10) s^cm" 3 . 

189 The differential multi-phonon production rate density is given by 

= n4 He <74 He P c — — s n (A)A. (11) 

190 ri4 Ho = is the number density and o"4 He = is the cross section per helium atom. E c and k c 

191 are the maximum kinetic energy and wavenumber of the neutrons trapped in the vessel. sji(A) 

192 is the scattering function sn(Q,u), evaluated on the dispersion curve of the free neutron, i.e. 

193 d = hQ 2 /(2m n ), with Q = 2vr/A. It is modeled^ to fit data from [28, 29] which is also shown 

194 in fig. 6a. Integration of dPu/dX from 0.52 nm to oo for the measured incident differential flux 

195 provides an estimated lower limit of 0.42 s cm - 3 for the integral multi-phonon production rate 

196 Ph. Difficult to estimate is the contribution to Pn from the wavelength range A < 0.52 nm. 

197 First, there is no direct experimental information available about sn(A) and second, this region 

198 has a strong weight due to the relatively large d(p/dX. The A-dependent measurements of UCN 

199 production by Baker and colleagues [T7J indicate that a multi-phonon contribution exists also 

200 from the region 0.4 nm < A < 0.52 nm. The measured integral multi-phonon and single-phonon 

201 production rates were found to add up evenly to the measured integral production rate |30j . 

202 This provides an experimental hint that for A < 0.4 nm the multi-phonon contribution might 

203 be negligible. The fit function shown in fig. 6 was therefore designed to extend the range of 

204 available data with a smooth decrease to zero at 0.4 nm. Integration of dPn/dA using the 

205 complete fitting function for sn(A) yields Pn = 0.56 s _1 cm -3 . Estimating the uncertainty of 

206 Pn as the difference between the values with and without inclusion of the contribution from the 

207 range 0.4 nm < A < 0.52 nm, we may expect a total UCN production rate of 



P th = (0.61 ± 0.10)i s^cm" 3 + (0.56 ± 0.14) n s^cm" 3 = (1.17 ± 0.17) s^cm^ 3 . (12) 

208 Evidently, more experimental information about sn(A) is required for a secure prediction of 

209 the UCN production rate induced by a white cold neutron beam. Note that this is particularly 

210 necessary for the region of short wavelengths, whereas the contribution to Pn from the region 

211 A > 0.8 nm is only 0.03 s -1 cm~ 3 . We also note that predictions from the two model calculations 

212 [20\ [2T] for the region A < 0.52 nm give mutually inconsistent results. 

2 More details about the modeling of su (A) will be described in a forthcoming paper. 
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Figure 6: a) Differential flux at the neutron guide NLf of FRMII (left scale), and measured 
scattering function sn(A) and fitted curve (right scale). The data point above the curve is due to the 
roton-maxon resonance in the multi-phonon dynamic structure factor, b) Calculated differential 
production rate, including single- and multi-phonon contributions, employing the measured d(/)/dX 
and the fit function for su (A) from fig. 6 a). 



213 5 Discussion and conclusions 

214 Our experiments demonstrate the feasibility of a versatile, intense UCN source at a cold neutron 

215 beam, providing UCN for experiments at room temperature. With the present prototype, at 

216 T = 0.82 K and using a small extraction hole with area A = 2 cm 2 , about 60 % of the UCN 

217 produced in the Fomblin grease-coated converter vessel reached the detector. Compared to our 

218 previous experiment with the uncoated electropolished stainless steel converter vessel [22], the 

219 count rates of accumulated UCN observed immediately after opening the UCN valve, and also 

220 the integral UCN counts were almost a factor three higher, despite the lower Fermi potential of 

221 the Fomblin wall coating. This increase of UCN output is due to longer storage time constants 

222 and a faster passage of the UCN through the extraction hole. Values for ta as derived from 

223 eq.© can be compared to the gas-kinetic equation 

4V 

ta = — A , (13) 
vA 

224 where we insert V ~ 2.4 1 for the volume of the converter vessel, and the mean velocity v of the 

225 trapped UCN, which, for a uniform distribution in velocity space, is given by v = 3t> max /4 with 

226 = Vp.waii — Vp 4jj e , where m n is the neutron mass. From the data for the Fomblin- 

227 coated vessel at T = 1.26 K we thus expect ta = 14.8 s. The measured value, ta = (27 ± 3) s, is 

228 much closer to the expected value than in our earlier experiments with the uncoated converter 

229 vessel. There, from eq. (!13p with the Fermi potential of stainless steel we may expect ta = 11-3 s 

230 but measured ta = (58 ± 13) s. A possible explanation is that UCN trajectories do not explore 

231 quickly enough the available phase space in the electropolished vessel, a hypothesis already 

232 raised earlier by W. Mampe and co-workers for their liquid- wall UCN bottle experiment |31j . 

233 Thus a certain roughness of the vessel walls seems to be necessary for efficient extraction of the 

234 UCN through a small hole. This should be taken into account in the design of a large converter 
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235 vessel in order to keep the time needed to extract the UCN reasonably short. 

236 The UCN production rate densities determined with two partly independent methods agree 

237 very well, eP\ = (0.80 ± 0.08) s^cm -3 , and eP 2 = (0.79 ± 0.09) s _1 cnr 3 . In our earlier 

238 experiment with the uncoated stainless steel vessel, we obtained the values eP\ = (0.91 ± 0.21) 

239 s _1 cm~ 3 , and eP 2 = (1.09 =L 0.25) s^cm -3 . The accuracy of the comparison has thus been 

240 improved by more than a factor two in the present experiments. More important, also the 

241 comparison with the theoretical expectation has been improved by virtue of a dedicated TOF 

242 analysis of the cold neutron beam spectrum, due to which we may expect the UCN production 

243 rate density Pth = (1.17±0.17) s _1 cm -3 , with single- and multi-phonon contributions of similar 

244 size. Since the detector and extraction efficiencies, expressed by the factor e, cannot be expected 

245 as perfect, the measured values come indeed very close to Pth- 

246 In the present experiments the saturated UCN density, normalised to the production volume 

247 V p did not exceed 46 UCN per cm 3 . To avoid activation of the stainless steel vessel only a quarter 

248 of the total volume was irradiated with the cold beam. Also for practical reasons, significant 

249 divergence losses due to beam collimation were accepted. In addition, for simplicity of these 

250 first tests, we did not yet employ a wall material with low neutron absorption and a high Fermi 

251 potential of 250 neV or beyond, such as Be, BeO or diamond-like carbon [32|, 133] . There one will 

252 gain in UCN density both due to the larger energy of the storable neutrons and due to a larger 

253 storage time. In order to produce a larger total number of UCN, the converter vessel can be 

254 made much longer, due to the low cross section of 27 x 10 -27 cm 2 for neutrons with wavelength 

255 0.89 nm [53], corresponding to a mean free path of about 17 m. With a properly designed 4 He 

256 converter and using an existing intense cold neutron beam at a high flux reactor, UCN densities 

257 in the order 10 4 per cm 3 with a total UCN number up to several 10 8 seem within reach. 
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